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Abstract 

It is argued that the Universe reheating in bouncing cosmologies could be explained via grav¬ 
itational particle production, as due to a sudden phase transition in the contracting regime. To 
this end, it is shown that, in the context of Loop Quantum Cosmology, gravitational production 
of massive particles conformally coupled with gravity in a matter-ekpyrotic bouncing Universe, 
where the sudden phase transition occurs in the contracting regime, yields a reheating temperature 
which is in good agreement with cosmological observations. 
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1 Introduction 

The issue of the Universe reheating in the matter bounee seenario (for descriptions of this scenario, see 
CD, via gravitational partiele produetion of light partieles minimally eoupled and nearly eonformally 
eoupled with gravity, has been reeently addressed, respeetively in |l2l and The idea is quite 
simple: to get effieient reheating one needs a non-adiabatie transition between two different phases, 
in order to obtain enough gravitational partiele ereation. This proeess is usually ealled preheating. 
At that time the Universe is far from thermal equilibrium, but then the created partieles decay into 
very light particles whieh interaet among themselves, thus producing a relativistic plasma in thermal 
equilibrium and evolving like radiation. Einally, when the energy density of this relativistie plasma 
starts to dominate that of the baekground, the Universe is reheated and enters in a radiation dominated 
phase which matches well with the eorresponding one for the hot Eriedmann Universe. 
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In the matter bounee seenario the non-adiabatie transition eould be produeed in the eontraeting 
phase. In faet, a transition from matter-domination (this phase is essential in order to guarantee seale 
invarianee of the power speetrum of perturbations, beeause modes that leave the Hubble radius during 
this phase have a flat power speetrum 01), to an ekpyrotie phase with equation of state P = up where 
u > 1 eould be assumed in the eontraeting regime. The model thus obtained is being ealled in the 
literature a matter-ekpyrotie bounee seenario [l5l, and sinee in the ekpyrotie phase the energy density 
of the field, namely p, evolves like whieh in the eontraeting phase inereases faster 

than a“®, anisotropies do beeome negligible. (Note that the energy density of the anisotropies grows 
in the eontraeting phase as a“®, that is, faster than the matter energy density, and thus, without an 
ekpyrotie transition the isotropy of the bounee would be destroyed; this is the so-ealled Belinsky- 
Khalatnikov-Lifshitz instability IH.) Moreover, the field energy density also grows faster than that of 
the partieles produeed, what means that the field dominates the Universe evolution in the eontraeting 
phase. But, when the Universe bounees, the energy density of the relativistie plasma generated by 
the deeay of the partieles ereated due to the phase transition will eventually dominate, and thus, the 
Universe will finally beeome radiation-dominated, in a way whieh may perfeetly mateh the results for 
the hot Friedmann Universe. 

This situation is similar to the one oeeurring in inflationary eosmology. To wit, for inflationary 
models with potentials not having a minimum (the so-ealled non-oseillatory models [|71), there is an 
abrupt transition from a quasi de Sitter regime to a kinetie-dominated period. During this transition, 
massless, nearly eonformally eoupled partieles are ereated whieh energy density evolves like pr ~ 
a“'^. On the other hand, when the inflaton field, <1), enters a kinetie-dominated period its energy density 
behaves like ~ [l8l|9l[l0l, what means that the inflation energy density deereases faster than 

that of the radiated partieles, and thus, the Universe will beeome radiation dominated, to eventually 
mateh the hot Friedmann Universe. On the other hand, in inflationary models with potentials whieh 
do exhibit a minimum, adiabatieity is broken when the inflaton field oseillates [fTTI ; at that moment 
massive partieles are ereated and, after the preheating stage, they deeay into light partieles whieh 
thermali / e and whose energy density eventually dominates the one of the inflaton. Therefore, the 
Universe beeomes radiation dominated, too. 

In the present analysis we will eonsider, in flat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
spaeetime, the matter-ekpyrotie bounee seenario in the eontext of Loop Quantum Cosmology ifT^ . 
where a phase transition to the ekpyrotie phase is produeed in the eontraeting regime. One of the 
most remarkable properties of this model is its simplieity, in the sense that in the non-singular bounee 
the adiabaeity is not broken, sinee it is well-known that given the Equation of State P = wp, dealing 

with Loop Quantum Cosmology, the seale faetor evolves as IfTSlI a(t) = + wYt^ -I-1 j , 

where pc is the energy density at bouneing time and Mpi the Planek mass. One ean easily see that, 
at bouneing time (t = 0), the seale faetor is an smooth funetion what preserves the adiabaeity. This 
situation does not happen—i.e., the transition from eontraetion to expansion is not adiabatie—in 
other bouneing eosmologies; for example, in eosmologies within the framework of General Relativity 
where the violation of the null energy eondition at bouneing time (needed in order to have a bounee in 
the flat FLRW geometry), eould be possible by ineorporating new forms of matter, sueh as phantom 
lfT4l or quintom fields ifTSll . Galileons [HHl, or phantom eondensates IITtI . 

In our simple study, we also deal with heavy massive partieles eonformally eoupled to gravity, be- 
eause in Fourier spaee, these partieles ean be depleted as a set of harmonie oseillators whose positive 
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frequency is time dependent. This allows to easily interpret the number of particles created and its 
energy density as the 00 component of the stress-energy tensor—where the vacuum zero point energy 
of the set of oscillators must be subtracted, in order to get a well-defined quantity (see, e.g., HU). 
Then, the energy density of the created particles will be calculated and, in the expanding phase, when 
it is of the same order as the background energy, the Universe will get reheated, with a temperature of 
the order of this energy density power 1/4. We will see that the reheating temperature will basically 
depend on the energy density of the background when the sudden phase transition is produced, and 
thus, by properly choosing this energy density, the reheating temperature will be shown to match well 
that of astronomical observations. 

Units used in the work are: h = c = 1. 


2 Particle production of massive conformally coupled particles 


Let 0 = f be a massive scalar field conformally coupled with gravity, a being the scale factor. In 
Fourier space, the corresponding Klein-Gordon equation, in the FLRW spacetime, is given by the 
following set of harmonic oscillators, with time dependent frequency [fT9l : 

x'k + (-jl{v)Xk = 0, ( 1 ) 


where ul{r]) = \k\‘^ + being rj the conformal time and m the mass of the field. 

Using the instantaneous diagonalization method COl . the square of the /9-Bogoliubov coefficient 
is given by 


mv)? 


1 


^ {Wkid)? + ^ki.d)\Xk{rii)?) 



( 2 ) 


We assume that at some given time, r]E, there is a sudden phase transition between two adiabatic 
phases. Note that, as ~ H^m (H being the Hubble parameter), if one assumes H < Hmax ^ Efi 
(i.e., the generated particles have a large mass), where Hmax is the maximum value of the Hubble 
parameter—what always happens in holonomy corrected Loop Quantum Cosmology because in the 
plane (FT, p) the modified Friedmann equation depicts an ellipse that is a bounded curve lf2T]| —then, 
for n = —1 + {r + l)a -f (s -f l)/9 with a + (3 — 1 > d, one gets 


dr/" \ J \ dp^ ) ’ 


in particular, for rr = 1, a = 2 and r = s = j3 = d, one has a;/ <C 

During the adiabatic regimes we use the first order WBK solution [l22l of ([I]) 


xiAd) 


, where lUi ^ 

V21Ui,,(r/) 


1 

/ l" 

3 

1- 

3 1 



4 

\0Jk) 


(4) 


together with its conjugate, to approximate the mode solutions. 

In the first phase, the mode corresponding to the vacuum state is approximated by xi,k- Before 
the sudden phase transition this mode becomes akXi,k + bkXi,k^ with 

bk = -i^[xi,k{VE)'^Xi,k{vt)], 


(5) 
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where W is the Wronskian and Xi,k(VE) (resp. Xi,fc(^E)) *^he value of the mode just before (after) 
the time tje- ^ simple ealeulation using the relation |afcp — = 1, yields 

l&(>))l^ = ^ [Ix'i.tWP + ^l{v)\xi.t{v)\‘‘] 

(a»6I.|(\''i,.(>)))" + i*'?(>))(Xix(>)))"l) = \h\^ + O (||||) = 16^1^ (« 

Assuming that the derivative of a has a discontinuity at some time value, p = Pe (namely that the 
Hubble parameter is discontinuous at r; = pe), one gets 


I**’')' =-- 

which, in terms of the cosmic time, reads 

\n /,\|2 ^ 

^ 16 (|A;p + m^a\Y ’ 


(7) 

( 8 ) 


where ue = a{tE) and 

When the Hubble parameter is continuous but its derivative is discontinuous, that is, the scalar 
curvature is discontinuous, at some instant fg, a similar calculation leads to 


\n /.\\2 ^ Re^ _ 

“ 2304(|fc|2 + m2a|)4 64(|fc|2 + m2a|)4 ’ 

where R{t) is the scalar curvature and = R{t%). 

On the other hand, the density of produced particles and their corresponding energy density are 
respectively given by [l23ll 


NM = 


\(3k{p)\^d^k] p^{p) = 


ujk{,p)\(dk{,p)\ d'^k. (10) 


(27ra(77))3 ’ {2Ea{p)fa{p) 

Note that, from the definition of the /9-Bogoliubov coefficients (O, we can see that the energy density 
pit) is the 00 component of the stress-energy tensor where, in order to obtain a convergent quantity, 
one needs to subtract the vacuum zero-point energy IIT^ . 

As a consequence, in the case when the Hubble parameter has a discontinuity at f = Ie, from 
formula ([8]) one obtains, in cosmic time. 


N^t) = 


miH+ - H 


■)2 

e) 


5127r 


Ue 

ait) 


Px(^) = 


m 


'■{Hj - H e? 

32ir2 


Ue 

ait) 


ait) 

O-E 


( 11 ) 


where 


h 





(a;2 -|- 1)3 


dx 



O-E 

a(0 _ ^ 

O-E 

CLE 


( 12 ) 
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From this result it follows that, when — 1, 

CLE 


Pxit) = 


- Hif (a(t) 


128 ir^ 


Qe 


and, consequently, if the sudden phase transition takes place in the contracting phase, then the gen¬ 
erated particles will evolve, in a part of the contracting and expanding phase, exactly like radiation 
does. However, if these particles do not decay as a relativistic plasma then, when , they will 

evolve as matter since, in this case 




/ ttE y 

5127r Va(t) / 


The same happens when the Hubble parameter is continuous but the scalar curvature has a dis¬ 
continuity at f In such case, from formula ®, one obtains, in cosmic time. 


Nxit) 


4096m7r \a{t)) 


Px(^) 


{Ht - He? 

1287r2 



(13) 


where /4 j 




(x2 + l)4 


dx, is related with 



through the formula 


h{y) 


1 d{y^h{y)) 

Qy dy 


In the asymptotic case ^ -C 1, one has 


Pxit) = 


{H^ - 
15367r2 



4 


and, when ^ > 1, 

' rj c' ' 


Pxit) = mN^it) 


PaEV 

40967r \CL{t) J 


Finally, it is important to stress that, when the Hubble parameter is continuous but the curvature is 
discontinuous, the energy density does not depend on the mass of the field. 


3 Bouncing models 


In this Section we will study two simple bouncing models coming from Loop Quantum Cosmology, 
where holonomy corrections modify Friedmann’s equation as follows lf24ll 




— (i - 
3M% V pj 


( 14 ) 
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being Mpi the reduced Planck mass and pc the so-called critical energy density, the value of the 
energy density at the bouncing time, its maximal value. From this equation it is clear that, when 
p pc, one recovers the standard Friedmann equation, that is, holonomy corrections can then be 
disregarded. Moreover, we can see that Eq. (fT4l) corresponds to an ellipse in the plane {H, p) llITI . 

As we have explained in the Introduction, holonomy corrected Loop Quantum Cosmology provides 
a non-singular bounce that preserves adiabacity. In fact, if during the transition from contraction to 
expansion the Universe is dominated by a barotropic fluid with Equation of State P = wp, the solution 
of (fT4l) will be given by ifTSll 

aW = (^^(1+ +1) , (15) 

which has an smooth behavior at the bouncing time t = 0, and thus, leads to an adiabatic transition. 

In both models we will assume that, in the contracting phase, there is a sudden transition from the 
matter domination phase to the ekpyrotic one, that breaks adiabacity. And since, as we have already 
seen, holonomy corrected Loop Quantum Cosmology preserves adiabacity at the bouncing time, in 
both models all particles are created, in the contracting regime, during this sudden transition. 


3.1 First model 


In this model we will consider a discontinuous Hubble parameter att = tp which scale factor is given 
by 


a{t) 



when t < tp, 
when t > Ie, 


(16) 


where we assume te > 1 in order to have an ekpyrotic phase. 

Note that for t < Ie, the scale factor is the solution of the standard Eriedmann equation H'^ = 
for a matter dominated Universe (P = 0), and for t > Ie, the scale factor is the solution of 
the holonomy corrected Eriedmann equation (1141) . for an ekpyrotic Universe whose EoS is P = wp. 
Assuming moreover that the sudden transition is produced at very early times, ^|^(1 + S> 1, 

one obtains and thus, from Eq. (fTTI) . the density of generated particles reads 




5127r(l -f wy \a{t) 


(17) 


Now we will show that this amount of produced particles is enough to reheat our Universe. To 
prove this statement note, first of all, that immediately after the transition, the energy density of the 
produced particles and the background energy density are, respectively. 


m^w'^p^ 


Pe 




( 18 ) 


where p^ is the background energy density at the end of the matter domination phase. 
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Then, in order to have a sub-dominant energy density for the partieles ereated at the beginning 
of the ekpyrotie phase—so that baek-reaetion effeets ean be disregarded—one has to assume that 

Mpi ^ raw. Sinee in the contraeting phase the baekground energy density increases as (i 


L PI I I HXJ . All lllW V-'V/lllA ClV-'llll^ J^llU-OW lllW L/ClV-'AV^l V/UllVl WllWA^^ VlWllOAl^J' AllV-'A WtlO WO CIO 1 (2(£) J 

it will be dominant during the whole phase, and only in the expanding one will both energy densities 
be of the same order, at some time t = tR. But then, the order of the reheating temperature may be 
obtained as Tr ~ p][^{tR). 

Secondly, as we have already explained in the Introduction, the particles created are far from 
being in thermal equilibrium. The more massives particles will decay into lighter particles, which 
will interact through multiple scattering, thus redistributing their energies until they yield a relativistic 
plasma in thermal equilibrium (see [|^ for a very detailed description of this process). 

Let r be the decay rate of massive x-particles. The decay will be accomplished at time when 
r|fdec — — 1- If, in order to simplify things (even though it is not essential), we assume this 

happens in the contracting phase when holonomy corrections can be neglected, then it turns out that. 


as for t > tE one has H = 


3(1-1-ui)t 


, one gets 




H, 


dec 


3(1 +in) H^H, 


rsj 


-1 


TJ r\j 

^dec — 


2Tm 


dec 


2T + 3(1 + w)H^ 


(19) 


This means that the decay rate must satisfy the constrain 


0 < \Hdec\ < 


Pc 


Mh 




( 20 ) 


because we are assuming that the decay is produced, in the contracting phase, whenever holonomy 
corrections can be disregarded. 

When the decay is effective, the background energy density reads 


Pdec = 


2TH^ 


2T + 3{1 + w)Hj 


Pa¬ 


ul) 


And thus, once the decay has taken place, a re-distribution of energies among the different particles 
occurs —kinetic equilibrium —and, also, an increase in the number of particles —chemical equilib¬ 
rium. That is to say, in such process, in order to obtain a relativistic plasma in thermal equilibrium, 
both number-conservating and number-violating reactions are definitely involved. 

Finally, to calculate the reheating temperature, one has to impose that both the energy densities of 
the relativistic plasma and of the background are of the same order. This will happen, in the expanding 
phase, when a{tR) 3> a^, namely 


p{^r) — 


Pe 


(1 + w )2 


a{tR)) Pxitn) + 


m'^w'^p^ 


ttE 

a{tR) 


( 22 ) 


what means that 


ttE 

a{tR) 


2 2 
m w 


rsj 


3w-l 


(23) 
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And thus, the reheating temperature will beeome 



(24) 


Note that in the asymptotie ease tu 3> 1, the reheating temperature has this simple expression 



(25) 


and in this ease, assuming from reeent observations that Tr ~ 10 (see [|3), one eoneludes that 


the energy density of the baekground at the end of the matter domination stage has to be neeessarily 
of the order 



That is, our model will definitely mateh observations, namely, the amount of created particles is 
enough to reheat the Universe, when the phase transition is produced at very low energy densities, or 
equivalently, at very early times in the contracting phase. 

To finish this example, it is instructive to compare it with its inflationary dual, that is, a Universe 


characterized by a sudden transition from the de Sitter phase to a radiation-dominated one. Suppose 


inflation is produced by a false vacuum energy density pp, then the scale factor is given by [j8l 



for T] < t]e < 0 
for T] > r]E, 


(26) 


where He is the value of the Hubble parameter during the de Sitter phase. Then, for heavy massive 
(m S> He) particles conformally coupled with gravity, since H{t]) is discontinuous at p = pp, one 
can use Eq. (fTTl) to obtain p^itE) ~ Thus, after the decay into light particles and thermal- 

ization, the reheating temperature, as in the matter-ekpyrotic bounce scenario with w 3> 1, is given 
by 



(27) 


which is another manifestation of the existing duality between the de Sitter regime in the expanding 


phase and matter-domination in the contracting one, pointed out, for the first time, in [|^ . 

3.2 Second model 

As second model of the matter-ekpyrotic bouncing scenario in Loop Quantum Cosmology, we will 
consider a sudden transition where the Hubble parameter is continuous but the scalar curvature has a 
discontinuity, at f = fp. If one considers Eq. (fT4l) with and EoS of the form 



0 when p < pE 
wp when p > pE-, 


( 28 ) 
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in the contracting phase, and P{p) = wp in the expanding one, the scale factor will be given by 


a(t) 



when t < tE, 
when t > ts, 


(29) 


where to = tE-^, being He = 


— ' 


the value of the Hubble parameter at the time of the phase 


transition. Note that, since we are assuming that the transition is produced at very early times, during 
matter domination, holonomy corrections could safely be disregarded, and thus, the exact value of the 

scale factor a{t) is approximately given by a^; (-HH. 


tE—to 


2/3 


In this case, R'p — = —9wH^, and then the density of produced particles is 


= 


4096m7r \a{t) 


Qe 


(30) 


On the other hand, at the transition time the energy density of these generated particles evolves as 


Pxit) ~ 


MU' 


(31) 


As a consequence, such energy density will be sub-dominant at the beginning of the ekpyrotic phase, 
provided w'^pe MU- And since the background energy density increases fastly in the contracting 
regime and also decreases fastly in the expanding one, the Universe will eventually reheat, in the 
expanding phase, when both these energy densities are of the same order. This will happen for 

O-itR) 3 > ttE- 

Imposing this condition, one obtains 


as ^ f w‘^Pe \ 

V mU ) 


(32) 


and thus, the reheating temperature will be 


T, 


R 


Pl/Utn) 



3(m+l)' 

Pe I 

fP eW^\ 

mu 

( MU ) 


1 1/4 


M. 


PI- 


(33) 


Finally, note that in the asymptotic case, tu 3> 1, the reheating temperature acquires the simple form 


Tr 


rsj 



1/2 

Mpi 


(34) 


and, in this case, assuming Tr ~ 10 ^Mp;, it turns out that, in order to match with recent observa¬ 
tional data, the energy density of the background at the phase transition must be of the order 


Pe ~ 


w 
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4 Some remarks on massless particle production 


When one deals with partieles non-eonformally eoupled with gravity, the diagonalization method 
results, in general, in a divergent number of produeed partieles. This means that, in sueh ease, the 
eoneept of partieles ereated has to be defined in a different way. We believe that, in the ease of 
massless partieles non-eonformally eoupled with gravity, the modes assoeiated to the vaeuum must 
be those whieh, at an early time, satisfy the asymptotie eondition Xk{v) 

Sinee the Klein-Gordon equation for massless partieles reads 


^-i\k\r] 


Xk + 




60- X. = 0, 

a 


(35) 


where ^ is the eoupling eonstant (equal to 1/6 when the eoupling is eonformal), with our definition, 
the modes that define the vaeuum state in the matter dominated phase are 


xr'\v) = e 


b'matt 


where H!)^2au 





(36) 


m 


is Hankel’s funetion Il27l . with Umatt = y | “ 12,^. Here, we will assume that 
order that Umatt beeomes a real number. Note that these modes are the same whieh define the vaeuum 


_3_ 

16’ 


in the de Sitter regime in the expanding phase beeause, in both eases, ^ this is the duality 


beetwen matter-domination during eontraetion and de Sitter phase during expansion 

On the other hand, if the phase transition takes plaee at very early times, when holonomy eor- 
reetions ean be safely disregarded, then the modes whieh define the vaeuum state in the ekpyrotie 
regime, just following the phase transition, are given by 


= e-'<- 


^ekpy 




(37) 


1 I 2(l-3m) 
4 (1-|-3«;)2 

eoeffieients are ealeulated as 


with Uekpy = 


(1 — 60- Consequently, the squared modulus of the beta Bogoliubov 


= 


w 




What is remarkable is that for te S> 1 one gets Uekpy = and thus, — 


p-i\k\v 


(38) 


That is. 


when re S> 1 the phase transition from the matter to the ekpyrotie phase in the eontraeting regime 
is exaetly the same as the phase transition from the de Sitter to a radiation dominated phase in the 
expanding one. This last phase transition has been studied in several works, in order to explain the 
reheating in inflationary models via gravitational partiele produetion (see, e.g., |[8l 1^12^11^ ). 

Then, for w 3> 1, when one eonsider modes well outside of the Hubble radius, i.e., whieh satisfy 
\kr]E\ "C 1, we ean use the formula for small arguments 


7/f (.) 9. - (./2rm - ^ (.nr 


( 39 ) 
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to obtain, after matching at rjE with the mode + /9fc-^7=, or using Eq. (l38l) . the formula 

^/ 2 \k\ \/ 2 \k\ 


— (|A^,e|/ 2)-"''-"“-' r2(^„„„) (1/2 - . 

lOTT 


(40) 


And in the opposite case, \krjE\ 3> 1—i.e., for models well inside the Hubble radius—from the 
asymptotic expression for large arguments [l27ll 





e 


—iz 


1 


(41) 


after matching airjE, one obtains 




(4^Lt/ - 1)^ 

m^E 


(42) 


It follows from these expressions that the energy density of the produced particles is always ultraviolet 
divergent, and also, that for iZmatt > 1 (i-e., ^ < ^) the number of particles created is infrared 
divergent, although their energy density is only infrared divergent for fmatt >3/2, i.e., for ^ < 0, 
what includes massless particles minimally coupled to gravity. 

To avoid infrared divergencies one must assume that the coupling constant ^ belongs in the range 
(5/48, 3/16) and, in order to have a finite energy density one needs to assume, as was shown in [[8]|, 
that the scalar curvature is regular (continuous, at the very least) during the phase transition. 

Simple formulas can be obtained if one only considers the production of particles whose modes 
leave the Hubble radius before the phase transition, i.e., for modes satisfying \kriE\ < 1. In this case, 
when 5/48 < ^ < 3/16 0 < Vmatt < 1 (to avoid infrared divergencies), using Eq. (l40l) one gets 


Nxit) = 




2567r3(l - 


r^(z/ma/t) 


2 ^matt 


2 


He 


3 



and 


Pxit) = 




2567r3(3 - 2VrnaU) 


■T {ymatt) I „ ^matt 



(43) 


(44) 


Note that the case Umatt = 1/2 ^ = | is the conformally coupled one, and thus, there is not 

particle production there. 

Einally, provided the phase transition is smooth enough, the energy density of the massless particles 
being produced will be of the order p^{t) ~ (| - ~ (I “ because 

for tu S> 1 the approximation j — 1 holds in the ekpyrotic case. On its turn, the reheating 

temperature will be of the order 


T, 


R 


2 ^matt 


Pe 

Wi 


1/2 


M, 


Ph 


(45) 


which coincides with the reheating temperature obtained in the transition from the de Sitter phase to a 
radiation-dominated one—because in this inflationary case one has the well-known result |l9l [TH [28]| 
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Px ~ (i “ (compare with Eq. (HU)), with i/ = Umatt = y/f — 12,^. This illustrates, onee 

again, the duality existing between the transition form the de Sitter regime to a radiation-dominated 
Universe, in the expanding phase, and the transition from matter domination to an ekpyrotie regime, 
in the eontraeting one. 


5 Conclusions 

We have shown in this paper that, for models arising in holonomy eorreeted Loop Quantum Cosmo¬ 
logy, a simple phase transition which takes place in the eontraeting phase, from matter domination to 
an ekpyrotie regime, ean lead in faet to the production of heavy partieles eonformally eoupled with 
gravity in a suffieient amount (see Eqs. (17) and (30)) and with enough energy density, able to reheat 
the Universe in the expanding phase. This is quite a remarkable conclusion. Moreover, as a bonus we 
see that the expression of the reheating temperature is aetually quite simple (see formulae (25), (34) 
and (45)), what does not happens in other bouneing models where adiabatieity is broken at bouneing 
time (see for instanee, Eq. (60) of flU, where produetion of minimally eoupled partieles and reheating 
are studied in the two-fields matter bounee seenario). In our results, the reheating temperature basic¬ 
ally depends on the energy density scale at which the phase transition oeeurs, what leaves us with just 
one degree of freedom in order to mateh this theoretical value with the observational data. This makes 
of the model a reasonably predietive one. Einally, we have emphasized that the reheating temperat¬ 
ure obtained in our model is the same one gets in inflationary eosmology via gravitational particle 
produetion, what elearly shows the dual relation that exists between the matter-ekpyrotie bouneing 
scenario in Loop Quantum Cosmology and the non-oscillatory inflationary models. 
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